Non-classical interference in integrated 
3D multiports 



Thomas Meany, 1 Michael Delanty, 1 Simon Gross, 1 Graham D. 
Marshall, 1 2 M. J. Steel, 1 and Michael J. Withford 1 

1 Centre for Ultra-high bandwidth Devices for Optical Systems (CUDOS), 
MQ Photonics Research Centre, Department of Physics and Astronomy, 
Macquarie University, North Ryde, 2109 NSW, Australia 
2 Centre for Quantum Photonics, 
H. H. Wills Physics Laboratory and Department of Electrical and Electronic Engineering, 
University of Bristol, Merchant Venturers Building, Woodland Road, Bristol, BS8 1 UB, UK 
thomas.meany@mq.edu.au 



Abstract: We demonstrate three and four input multiports in a three 
dimensional glass platform, fabricated using the femtosecond laser direct- 
write technique. Hong-Ou-Mandel (HOM) interference is observed and a 
full quantum characterisation is performed, obtaining two photon correlation 
matrices for all combinations of input and output ports. For the three-port 
case, the quantum visibilities are accurately predicted solely from measure- 
ment of the classical coupling ratios. 
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1. Introduction 

Optical multiports provide the ability to produce an arbitrary unitary transformation of a set 
of optical modes Multiports are therefore an extremely important component in optical 
quantum information processing (QIP) and state preparation 013), where producing entangled 
states is a key motivator for such devices flUO. Multiports and their applications in quantum 
optics were originally explored using bulk optic beamsplitters and phase shifters @, for which 
non-classical interference was observed in both three- and four-port splitters. Following from 
this work a study of a three-port fiber coupler was undertaken [7]. However neither of these 
approaches provide the scalability and stability requirements needed for large scale optical 
quantum circuits. The extreme difficulty of working with large numbers bulk optical elements 
and fibre components has limited the development of these devices. 



The challenges of scaling and stability in various quantum optical applications have been 
the main drivers for the recent explosion of interest in integrated quantum devices, such as inte- 
grated waveguide circuits [8 9 1, waveguide array quantum walks |fT0l[TTl[T2l and single photon 
sources lfT3l[T4l . Desirable qualities such as reconfigurability have been demonstrated [15], no- 
tably in lithium niobate-based devices which allow for telecommunications wavelength photon 
manipulation at GHz speeds [16|. Miniaturisation using silicon-on-insulator nanowire based 
photonic circuits have been shown ifTTIl and there is also great promise to incorporate these 
platforms with existing on chip detectors 1 18 1. In addition the integration of photon pair sources 
and circuits have been achieved using a nonlinear waveguide array fabricated in lithium nio- 
bate 11191 . The advantages of integrated components have been applied to the construction of 
integrated quantum multiports. Peruzzo et al. have demonstrated multiport functionality with 
an on-chip 4-port multimode interference device [20 1 in the silica-on-silicon platform. 

The majority of these demonstrations have relied on mature planar platforms, but there are 
limitations associated with the inherently two dimensional (2D) nature of planar quantum cir- 
cuits. For instance in the realm of optical quantum simulation |2TI . there are interesting Hamil- 
tonians which can only be effectively simulated using waveguides in three dimensional (3D) 
arrangements El . Moreover, as circuits become more complicated, both the circuit area and 
the proportion of "neutral" elements, whose only purpose is to provide transparent waveguide 
crossings, will become significant. Hence techniques which allow the fabrication of structures 
in three dimensions have some specific and unique possibilities. Such a technique is the fem- 
tosecond laser direct-write technique (FLDW) which allows for fabrication of waveguides in 
glass subsrates in 3D. This technique has been used to produce a number of devices show- 
ing quantum operation including directional couplers [23 1 (evanescently-coupled waveguides 
which act as beamsplitters), a 6-port continuous quantum walk on an elliptical array l24l . a 
discrete quantum walk on a line using entangled photons [25 1 and investigations of quantum 
interference in 2D arrays [26|. This has demonstrated the potential of the FLDW technique to 
fabricate useful structures for quantum applications. 

In this paper we implement a 3-port and 4-port device using a 3D FLDW technique. Each 
device highlights a different advantage of the 3D geometry. The 3-port is implemented as a 
single coupling region that replaces three 2-port devices JT). The 4-port is composed of four 
2-port directional couplers, but avoids a crossing that would require an additional coupler in 
2D. Taken together these kinds of structures could reduce the complexity and sensitivity of 
higher order concatenated multiport based circuits. This fact is established by our observation 
of 2 photon interference in both structures. 

2. Implementing multiports 

The fundamental idealised multiport device is the familiar beamsplitter which implements a 
desired 2x2 unitary operation. This can be generalised to a multiport operation on N input 
modes into N output modes L>n described by the unitary N xN matrix [6| 
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The integrated photonics analog of the basic 2-port beamsplitter is the directional coupler, in 
which two waveguides are brought into close proximity so that their evanescent fields can 
overlap and exchange energy [8 27 1. The degree of coupling is determined by the waveguide 
separation and interaction length which can be adjusted to produce directional couplers of any 



reflectivity. Using the well-known decomposition of the N x N unitary into a concatenation of 
beamsplitters and phase shifters 0"), it is possible to construct an arbitrary multiport in two 
dimensions in integrated optics by designing a suitable arrangement of directional couplers and 
phase delays. However as the number of directional couplers scales with A^ 2 , at large N, such a 
device is prohibitive to build due to fabrication errors. 

Previous investigations have looked at cases of bulk optic beamsplitters [6] , fibre optic di- 
rectional couplers Q, and multimode interference devices [20 1, to produce multiports In the 
next section we apply these techniques to the specific cases of a three and four port beamsplitter 
based on an arrangement of coupled waveguides 

2.1. Mode evolution in the three port ( fritter) 

In three dimensions a 3 x 3 multiport or tritter can be constructed from a single coupling el- 



ement involving three waveguides (see Fig. la I. This single element can replace up to three 
directional couplers and phase shifters [1|. The coupling coefficient G(z) describes coupling 
between waveguides 1 and 3, and g(z) describes the coupling between waveguides 1 and 2 and 
waveguides 2 and 3. In general, all three of the couplings could have different values, but, as 
will be discussed later, the 'isosceles' geometry considered here is appropriate as a result of the 
laser writing process we used. 



g(z) g(z) 

!# G(z) l 3 





(a) 



(b) 



(c) 



Fig. 1 : (a) Idealised coupling region of a tritter device. Writing considerations induce an effec- 
tive isosceles geometry with two coupling strengths g(z) and G(z). (b) A schematic of a 3D 
implementation of a tritter, where three waveguides taper from a planar arrangement to a trian- 
gular interaction region and taper back to a planar array and (c) the layout of the four directional 
couplers and the waveguide avoided crossing in the 4-port device 



Such a structure is described by a Hamiltonian of the form 
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The propagation constant j3 = (o/v p of the three waveguide modes is assumed to be common 
(v„ is the phase velocity and a>/2n the optical frequency). Explicitly including the common 
phase velocity v p in H means that operator evolution with distance z rather than time t follows 
naturally. 

After propagation through the interaction region, photons have encountered the mode trans- 
formation 



b = f/a 



(4) 



where the input and output mode vectors are a = (<Zi , a% , &■$ ) and b = (b\ , bz, 03 ), and the trans- 
fer matrix U is the solution of Heisenberg's equation, da/dz = — iC(z)a. The transfer matrix U 
completely determines the operation of the device, in particular the classical output intensities 
and quantum interference between photons in the device. 

2.2. Mode evolution in the four-port 



In contrast to the single element tritter the four port (depicted in Fig. lc 1 is made up of four 
directional couplers. Here we have used the 3D advantage to cross two waveguides without 
interaction, reducing the number of required directional couplers by one. As the device is made 
up of four discrete elements, the device can be modelled as a product of beam splitter operators. 
Assuming each directional coupler is identical, the operator corresponding to the four port 
transfer matrix is, 



U = ^2,4(l?)^l,3(r/)A(^)^(^)B3,4(T/)^l, 2 (T?) 

where the directional couplers have the beam splitter action, 
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and Pj(<j>) describes the phase difference between the inner crossing arms and the outer 
"straight-through" arms (See Fig. lei where Pj((j))dj = e^'^dj. Using these relations we find 
the transfer matrix corresponding to the four port operator Q is 
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3. Fabrication methods and design 



The devices were fabricated using the femtosecond laser direct-write (FLDW) technique which 
employs a tightly focused femtosecond laser to inscribe a localized refractive index change 
in glass. By translating the glass sample in the (x, y, z) directions with respect to the incident 
laser, arbitrary 3D regions of net-positive refractive index change (which act as waveguides) 
can be produced inside the host glass. The laser used to fabricate the structures is a titanium 
sapphire oscillator (Femtolasers GmbH, FEMTOSOURCE XL 500, 800 nm centre wavelength, 
< 50 fs pulse duration) with 5.1 MHz repetition rate. A 100X oil immersion objective (Zeiss 
N-Achroplan, NA = 1.25) was used to focus the laser inside boro-aluminosilicate (Eagle 2000). 
The sample translation was completed using Aerotech high precision motion control stages 
with 10 nm precision at writing speeds of 2000 mm/min. The waveguides were written us- 
ing pulse energies of 28 fii which created waveguides suitable for single mode operation at 
800 nm. The combination of low pulse energy, high translation speed and high repetition rate 
created waveguides which were at the at the very onset of cumulative heating modification 
mechanism [28 1 (implying that pulses arrive before the characteristic thermal diffusion time of 
the material). This results in smooth refractive index change and short device production time, 
allowing the fabrication of complex photonic circuits. The refractive index change is ellipsoidal 
in nature and it has been shown that laser written waveguides display non-isotropic coupling 
due to waveguide asymmetry [29|. Since the tritter is a triangular arrangement of modes (as 



shown in Fig. la 1 there are two effective coupling ratios, g(z) 7^ G(z). In both cases the waveg- 



uides are initially separated, in a planar array, by 127 fim before tapering an interaction region, 



where evanescent coupling occurs with typical spacings of 15 /im and subsequently tapering 
out to 127 /im planar spacing. This is best illustrated in Fig. [T[ which shows the waveguides 
tapering towards interaction regions and returning to the original spacign at the output of the 
device. We set out to fabricate a tritter with equal splitting between output ports of 33/33/33 
and a 4-port with equal splitting between output ports of 25/25/25/25. Hence we fabricated a 
range of devices with different length interaction regions targeting the ideal coupling ratios. We 
performed a classical characterisation to determine which structures best matched our intended 
constraints. 



4. Classical characterisation 



In order to determine the coupling ratios in the tritter and splitting ratios for the four port a 
classical characterisation was performed. If classical light of intensity M is injected into port /', 
the output intensity N at port k can be modelled as, 



N k = gj'epMlt/,--' 2 



where the input and output losses at port / are characterized by the pre-factors, < £ ; 
By forming the following ratios, 
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it is possible to relate the measured classical intensities directly to the unknown tritter couplings, 
g(z) and G(z) or four-port splitting ratios and cancel the loss terms. 

For the tritter we fitted a model to the classical measurements by assuming uniform coupling 
over an interaction region of length L, i.e. g(0 < z < L)L/v p — g and G(0 < z < L)L/v p = G 
with zero coupling elsewhere. Using maximum likelihood estimation [30], we varied g and G 

1 2 

in the transfer matrix U to find the best fit to the three experimentally determined fractions F l ' , 

13 2 3 

Fj 3 and F 2 3 . This involved minimizing the difference between the experimentally measured 

§12 13 2 3 

for Fj 2 , -Fj 3 and F 2 3 , weighted by the 

uncertainties in the fractions (found from the uncertainties of the measured intensities Nj) [?]. 

This fit determined g « 0.81 V p /L and G ~Q.5\v p /L, indicating that waveguides one and two 

couple more strongly than waveguides one and three. Using this fit we obtained the matrix of 

classical output intensities 



0.37 0.41 0.23 
0.41 0.19 0.41 
0.23 0.41 0.37 



(10) 



The ideal symmetric tritter has |C/,- ,-| 2 = 1/3, therefore we predict our device has a slightly 
asymmetric power splitting (6) . 

Similar to the tritter, we can use a maximum likelihood abroach to fit tj, the beamsplitter 
reflectivity, from tl i using the following six fractions, F 1 ' 2 , 3 , Fj \ , F 2 3 , F 2 \ , and F 3 \ . This 
fit determined r\ =~0.377. Substituting this fit into Eq. |7]i we find the matrix of classical output 
intensities for the four-port 
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Our four-port device differs from the ideal symmetric four-port, \Ujj\ 2 = 0.25 and has a non- 
uniform power splitting |6|.This is due to the fact that the value we measure for 77 = 0.377 is 
not the ideal value of tj = 0.5. 



5. Quantum characterisation 

We performed two photon characterisation on the devices using 804 nm photons produced 
using type-I spontaneous parametric down conversion (SPDC) source, as described in ll23l . 
The setup used a Toptica 402 nm blue diode laser which was focused into a 1 mm thick BiBO 
crystal cut to give a 6° opening cone angle. Down converted photons were produced at 804 nm 
in an output cone and then coupled into polarisation maintaining single mode fibres. These were 
then butt-coupled to the device under test using 127 jim spaced V-groove arrays. The outputs 
were then monitored using silicon avalanche photodiode detectors (SAPD) from Perkin Elmer. 
The 2-fold coincidences across all the output modes were measured simultaneously (in a 5 ns 
window) using a 4 channel time tagging unit. The arrival time of the photons into the chips was 
varied as a means of continuously controlling the degree of indistinguishability between the two 
photons and thereby moving between classical propagation of light and quantum interference. 

We detected signatures of two photon interference in the devices through the Hong-Ou- 
Mandel (HOM) effect 0T1 . whereby the probability of two photon coincidences are reduced 
(HOM dip) or enhanced (HOM peak) due to quantum interference in the device. The degree of 
quantum interference is quantified using the visibility, Vjj, defined for injecting a single photon 
into waveguides i,j and detecting in waveguides k,l as, 

where, the quantum and classical coincidence probabilities are, 

Qn = T^rWiMUji+UaUj,^ 2 (13) 

C-j = lU^Ujjf + lUtjUjtl 2 . (14) 

We see that the quantum two photon coincidence probability Q.'- is strongly dependent on the 
relative phases of the elements of U . When these phases destructively interfere, two photon 
coincidences are reduced compared to the classical case, leading to a large positive visibility (a 
HOM dip). Whereas when these phases constructively interfere, two photon coincidences are 
enhanced compared to the classical case, leading to a large negative visibility (a HOM peak). 

5.1. Visibilities in fritter 

As the visibilities are strongly dependent on the relative phases of the matrix elements of U, the 
tritter visibilities depend on the length and geometry of the coupling region as defined by the 
coupling functions g(z) and G(z) in Eq. 0. To illustrate the variety of tritter visibilities that can 
be observed in our 3D device, we plot the visibilities as a function of effective coupling gL in 
(Pig. |2j assuming uniform coupling over the interaction region. Figure [2a] shows a symmetric 
tritter as gL is varied. Only at one coupling ( the symmetric 33/33/33 case) do we see equal 
visibilities of 50 
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Fig. 2: Tritter visibilities as a function of effective coupling, gL. We assume uniform coupling 
over an interaction region of length L, g(Q < z < L)L/v p = g and G(0 < z < L)L/v p = G with 
zero coupling elsewhere, a) symmetric coupling region (G = g). b) Isosceles coupling region 
(G = 0.6234g). Vertical line corresponds to the fit to our experimental tritter with gL « 0.81 V p . 
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Fig. 3: Tritter HOM dips for (a) input ports 1 and 2, (b) for input ports 2 and 3 (c) for input ports 
1 and 3. (d) Comparison of measured tritter visibilities (error bars) and theoretical predictions 
(columns) from the simple tritter model in section |4| 



We measured nine HOM dips for the tritter (Fig. [3]) by injecting two photons into the device 
at different delays and observing coincidences at the output ports. The measured visibilities 
(Fig. 3d I were found from the relative depth of each dip. We see a clear reduction in the visi- 
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In section|4]we fitted a simple model of the tritter couplings g(z) and G(z) using only classi- 
cal intensity measurements of the device. As this fit uniquely determines U we can predict the 



visibilities for this model using Eq. ( 12 1. Therefore we are able to predict, using only classical 
intensity measurements measurements, the quantum correlations for the device. This is a conve- 
nient tool, since typically in order to produce a device with the correct coupling ratios a range 
of device must be fabricated around this indended parameter. It is much more convenient to 
perform a range of intensity measurements prior to any quantum measurements to confirm the 
device has the desired functioning. The predicted visibilities are compared to the experimental 



visibilities in Fig. 3d where we see a qualitative match, issues which are causing the minor 



mismatch between this prediction and the experimentally observed values are the assumptions 
uppon which the intensity dependent measurement is based. The model assumes a constant in- 
put intensity which can fluctuate during the course of the measurement and input output facet 
losses which may change also. Furthermore the model assumes the interaction region is con- 
stant in z but input and out tapering regions do contribute to the overall coupling which was not 
accounted for. In addition since the maximum likelihood measurment is a numerical simulation 
it has it's own inherent errors I30ll32l . 

5.2. 4-port 

We also measured HOM interference for all input and output combinations in the 4-port device. 



These measurements are shown in Fig. 4a 4f| We see a variety of different magnitudes and 
signs of the visibilities for the various input/output combinations. This is due to the quantum 



12 1, which exhibits constructive or destructive interference, depending on the 
me elements of U. The fact that all outputs had a visibility greater than 50% 



term Q^j in Eq. ( 
relative phases o 

(or less than (—50%)) indicates that the output states are strongly non-classical. 

We now compare the experimental visibilities to those predicted from the 4-port model in 
section [4] In contrast to the tritter the classical intensity measurements in the 4-port are insuffi- 
cient to determine the phase (j) due to the additional path lengths of waveguides two and three. 
Therefore after fitting 7] = 0.377 from the classical intensity measurements we used visibility 
measurements from a single input, namely inputs two and three, to fit (p. Using maximum like- 



lihood estimation we fitted <p from the six HOM dips in Fig. 4d This fit determined = O.Oln. 



After substituting the determined 77 and into the transfer matrix Eq. (|7]l we can use ( 12 1 to 
predict the other 30 visibilities in Fig.|4] The predicted visibilities are compared to the measured 
visibilities (see Fig. 4g 1 and display a qualitative match, correctly identifying the sign and mag- 
nitude of each visibility. We note that using other input combinations to determine (j) produced 
similar predictions. 
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Fig. 4: Four-port HOM dips for (a) input ports 1 and 2, (b) for input ports 1 and 3 (c) for input 
ports 1 and 4 (c) for input ports 2 and 3 (d) for input ports 2 and 4 (e) for input ports 3 and 
4. (f) Experimental four-port visibility matrix determined from the HOM dips, (g) Comparison 
of measured four port visibilities (error bars) and theoretical predictions (columns) from the 

frvnr_r\rvrt mr»H^l 



6. Conclusion 



Two photon correlation measurements were performed on integrated 3 and 4 port devices. 
The visibilities observed for our tritter compare extremely well with bulk optic demonstrations 

12 13 23 

where visibilities, for only one pair of input ports, Vj 2 ~ 50%, Vj 2 ~ 48% and V l 2 ~ 26% are 
observed (6|. In addition our tritter displays greater symmetry than a fibre optic based example 

12 13 

where visibilities, observed only for a single pair of input ports V x \ ~ 50%,^ , ~ 30% and 

Vj 2 2 ~ 25% 0. Our four port has shown a uniform output very similar to a bulk optic example, 

where visibihties of y/f - 75%,V, 1 2 3 - 66%, v/:, 4 66%, - -80%, ~ 75% and 

Vj 3 2 4 ~ 66% were observed |6] which out performs a multimode interference device described 
in 1 20 1 . This demonstrates the applicability of laser written photonic circuits for some quan- 
tum logic applications where a 3D capability is advantageous for circuit size or complexity 
requirements. 
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